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Abstract 


In this paper, the general problem of a shell containing a 
through crack in one of the principal planes of curvature and 
under general skew- symmetric loading is considered. By employ- 
ing a Reissner type shell theory which takes into account the 
effect of transverse shear strains, all boundary conditions on 
the crack surfaces are satisfied separately. Consequently, 
unlike those obtained from the classical shell theory, the angu- 
lar distributions of the stress components around the crack tips 
are shown to be identical to the distributions obtained from the 
plane and anti-plane elasticity solutions. Extensive results are 
given for axially and circumferentially cracked cylindrical shells, 
spherical shells, and toroidal shells under uniform in-plane 
shearing, out of plane shearing, and torsion. Taking advantage 
of the fact that the problem is formulated for "specially" ortho- 
tropic materials, the effect of orthotropy on the results is 
also studied in some detail. 

1. Introduction 

In recent years, a great deal of effort has been devoted to 
the study of cracked plates and shells. The reason for this 
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stems from the fact that a great variety of advanced structures 
such as aerospace vehicles, pipelines, different components of 
nuclear reactors, etc. are designed from the viewpoint of frac- 
ture mechanics. The early studies of cracked shells employed 
the so-called "classical" shell theory. This is an eighth order 
theory, and consequently can accoiiinodate only four boundary con- 
ditions on each crack surface. To make the number of unknown 
functions arising from the solution of the differential equations 
compatible with the number of independent boundary conditions, 
the transverse shear and the twisting moment are combined as the 
"effective transverse shear" which in turn is used to satisfy the 
boundary condition regarding the transverse shear, (See, for 
example [1-7], and for review and references [8]). The short- 
comings of the classical shell theory are well-known from the 
plate and shell solutions. The angular distributions of the 
moments around the crack tips do not conform to those obtained 
from the plane elasticity solution and the transverse shear has 
a strong singularity of order of -3/2 which is not physically 
acceptable. It has been shown that in plates, if one uses a 
Reissner type plate theory [9-12] which adequately incorporates 
the effect of transverse shear strains, then the angular dis- 
tributions become identical to those obtained from the plane 
and anti-plane elasticity solutions [13-18]. Motivated by the 
success achieved in solving the plate problems, in recent years 
a Reissner type shell theory [19-20] has also been applied to 
cracked shells [21-26]. However, in all these studies, the 
loading- is symmetric (i.e., the cracked shells are either 
under tension or bending). In practice skew- symmetric 
loading is also very important. Considering also the fact 
that the existing solutions are limited, the need for a 
refined solution and additional reliable results does not 
need elaboration. In this paper the general problem of cracked. 
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thin and shallow shells under skew-symnietr1c loading will be 
formulated by using a Reissner type shell theory. After 
analyzing the asymptotic behavior of the stress resultants 
around the crack tips, extensive results for circuit .rentially 
and axially cracked cylinders, and spherical and toroidal shells 
under general skew-symmetric loading will be given. The effect 
of material orthotropy on the results will also be studied. 

2. Formulation of the Problem 

The problem under consideration is that of a shell of 
thickness h containing a crack of length 2a in one of the princi- 
pal planes of curvature and subjected to skew-symmetric loading. 
The problem is formulated under the following assumptions: a) the 
shell is shallow and thin, b) the effect of transverse shear 
strains is included, and c) the material is "specially” ortho- 
tropic, i.e. the elastic constants are related as follows; 

where 


E = , V = 

Condition (2.1) will enable us to factorize the differential 
equations. Referring to figure 1 for the geometry of the shell 
and notdCion, the equilibrium equations can be written as: 


ij.j ' ° 

(2.2) 

i.i +q(X,.X2) = 0 

(2.3) 

'ij.j - Vi ' 0 (i’j = ’’2) 

(2.4) 
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where the indicia! notation and the summation convention are 
used, and (i,j«l,2) denote the stress, 

moment and transverse shear resultants, respectively, q(X-;,X 2 ) 
is the load normal to the X 1 -X 2 plane, and Z(Xi»X 2 ) is the equa- 
tion of the shell. Let Uj, W and 6^* (i“l,2) be respectively the 
displacement components and the angles of rotation of the normal 
to the shell surface. Using the generalized Hook's law, in terms 
of the displacement derivatives, the strains are obtained as: 




(2.5) 


Defining a stress function F by 
^ij “ ®ik ®jl ^kl 

e||^ being the permuation symbol, the equilibrium eqs. (2.2) 
are satisfied identically. Eliminating U-j and and substituting 
(2.6) into (2.5), we obtain: 

®im®jn®kp®lq®ijkl*^*mnpq ^ ^*ij®ik®jl^'**kl “ ° 

For "specially" orthotropic materials satisfying condition (2.1), 
the strain-stress resultant relations may be expressed as: 

^11 " W (^n/c^-vNgg), e-|2 TT “^IZ’ "^22 hf ^^^'^22"''^11 ^ 

( 2 . 8 ) 

^11 ‘^(^^®l,l'^^2,2^‘ ^12 "" 
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cu3 

^22 - ) . D - 120-v^T ’ 

(.2.9) 

cW “,1 h ’ W Vg - W.g + $2 > B - izn+v) 

(2.10) 

where c = (E-i/Eg)^. 

Using the relations (2.8)-(2.10) , the normalized quantities 
defined in Appendix A, the curvatures defined by 

a^Z 1 3^Z _ 1 a^Z _ 1 

R^j , Rg * ax^aXg R^g * 

(2.11) 

and the following new functions, 

93. 

n(x,y) = -jy 3^ . 

(2.12) 

as as 

tp(x,y) = + -^) " w , 

(2.13) 

equations (2.7), (2.3) and (2.4) become: 


- JZ -jp- - 2 ^)W - 0 

+ x2(1-kv2)(x^2 3^ - 2X^2^ sxsy ’‘2^ alP'' 

(2.14) 

= X'^Cl-icV^) ^ q 

(2.15) 

ic v^jjj - Ip - w - 0 

(2.16) 

~ n = 0 

(2.17) 
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The problem is thus reduced to the solution of the differential 
equations (2.14) - (2.17) under the boundary conditions as yet 
to be specified. 

For the crack problem the transverse load q*0, and noting 
that the crack is along a plane of principal curvature the 
shell parameter x -|2 ® 0. Further it is assumed that through a 
proper superposition, the problem is reduced to a perturbation 
problem where the only external loads are those acting on the 
crack surfaces. 

Defining 

eqs. (2.14) and (2.15) yield: 

+ 7? vf (1 -k72)<|) = 0. (2.19) 

A A 

Now assuming <j) in the following form, 

00 

<l>(Tt,y) = ^ g(x,a)e"''^“ da (2.20) 


and using the regularity condition at x=+«, (2.19) gives: 

4 m.x 

Z Rj(a)e , X > 0 

g(x,a) = < (2.21) 

8 m.x 

Z R4(a)e , X < 0 
. 5 <J 


Rj, (d=l,...,8) are unknown functions and m-j,...,m0 are the 
roots of the following characteristic equation: 
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m 2 a p + c (2 

- K ^2^ ■** (2kXi^X2^“^ ** 2K:X2‘^a^ + \2^)P^ 

+ (2kX-|^X 2^“^ “ KX2‘*a^ - + 2X^ - 2X.|^X2^)a^p 


+ (Xg^ - X^2)2 - 0 

(2.22) 

such that. 


RG(mj) < 0 , ~ J ” 1 » • • » »4 

(2.23) 

Similarly assuming, 


*1 * 

w(x,y) = ^ f(x,a)e'“^^“ da 

-00 

(2.24) 

03 

J2(x,y) = [ h(x,a)e“^'^“ da 

.... j 

—00 

(2.25) 

-00 

»|^(x,y) = ^ e(x,a)e"'’^“ da 

0 

(2.26) 


■■00 
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h(x,a) 


and 


J A^(a 

I AgCa) e 


r,x 

) e , X > 0 


-r^x 


, X < 0 


(2.28) 


e(x,a) = < 


R,-(a)p^2 m-x 

e J , X > 0 




(2.29) 


o 8 Ri(a)Pi-4 

X2 2 7 ■ n rV7 ' V"Zm\ V 2 2T 6 , X < 0 

5 U'Pj-UUg nij “ J 


where 


. - /■ 


a 


2 + 


■“Tt r 

K\ i-v; 


(2.30) 


Considering the symmetry of loading and geometry with respect 
to the y-axis, the stress resultants satisfy the following sym 
metry conditions: 


=* -Nj^j^(“X,y) , N^y(.x,y) 
Mj^j^(x,y) - -M^j^(-x,y) , Mj^y(x,y) 


Nj^(-x,y) 

Mj^y(-x,y) 


Vj^(x,y) = Vj^(-x,y) 


(2.31) 


In solving the problem, it is therefore sufficient to consider 
the x>0 portion of the shell only. The boundary conditions of 
the problem can then be expressed as follows: 
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Nj^j((0,y) = 0 

-o> < y < 

GO 



M^^(0,y) = 0 

-CO < y < 

CO 



lim 

x->0-^ 

Nxy(x,y) 

r— 

ti- 

ll 

-✓c < 

: y < 

✓c 1 

> 


v(0,y) = 

0 

ly| 

> /o 

J 

lim^ 

X-.0+ 

Mj^y(x,y) 

= Fg(y) 

-/c < 

■ y < 

> 


3y(0,y) ' 

= 0 

|y| 

> )/c 


lim 

x-^"^ 

Vj^(x,y) = 

= F3(y) 

-^c < 

: y < 

^ 1 


w(0,y) = 

0 

|y| 

> c 

j 


(2.32) 

(2.33) 


(2.34) 


(2.35) 


(2.36) 


where F-j(y), p 2 (y), p 2 (y) define the loading on the crack sur- 
faces and are known functions. Using eqs. (2.6), (2.8)-(2.10) , 
(2.12), (2.13) and (2.20)-(2.30) the stress resultants may 
be expressed as: 


Nxx(x>y) 


2tt 


4 m»x 4 w- 
E R.(a)e e da 
1 J 




JL 

2ir 


CO 

f 4 


•t ni.x 

i m^. R.-(a)e ^ da 

1 j J 


Nj^(x,y) 


1 

2^ 


4 rn.x _■ 

a E m.| R..(a)e ^ da 

1 j j 


(2.37) 


(2.38) 


(2.39) 
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a 1 


f” ^ 4 Rj(a)p^2(m^2_va2) 




,V e-iy« 


da 




TiX 


ar-jA^(a)e ^ da 


(2.40) 


Hyy(x.y) = 


4 R4(a)p,^{vlTlj2-a2) 


^^3^ e-^y« 


da 


i a (1-v)^ 


a r.jA^(a)e ^ e“^y“ da 


(2.41) 


R^(a)P^^ni^ e"’^V^^“da 


Mj^(x,y) = - ^ hxZ* J 1 (Xg^m^x^ViiKpj.-n 


. J.^ ^ ilplii f (c.2+r,2)Ai(o)e''’’' e‘^^“ do 


(2.42) 


iiir HA' 


j 

-00 


m.x 


. {<” 4 Rd(ot)P--^ni. ■•..■«> 

Vx(x,y) “ ^ j 5 ^ ^ 


- i ^ I a A^(a) e ^ da 


(2.43) 


Vy(x.y) - - ^ 


e-y« da 


“ f ( Xg^m . ^ ” X ^ ^a=^ ) U P j “ • y 


- K 


1-v 


*^1^ “iya 


A-j(a)e e 


da 


(2.44) 


Application of boundary conditions (2.32)-(2.36) would lead to a 
system dual integral equations. However, if one defines a set 
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Of new functions 1n terms of displacement and rotation deriva- 
tives, the problem can also be reduced to the solution of three 
simultaneous singular integral equations. 


Let, 




G-j(y) * 

» lim, 

av ^ 

3y ’ ^ ^ 

(2,45) 

G2(y) ' 

= lim 

3y 

(2.46) 

G3(y) ’ 

= lim, 
x-tO+ 

3W 

3y 

(2.47) 


Using the basic equations resulting from the formulation of the 
problem, the new unknown functions may be expressed as: 


Gi(y)=A f ? da 

I b. tl I I U 


(2.48) 


; I yj M 

•mca 


T r 4 Ri(ct)p..^ “iya ^ 


X 1 1-v 

2tt 2 


a r^ A-j(a)e da 


(2.49) 


G3(y) = - 


2ir 


4 Ri(a)Pi^ -iVa ^ 


(2, GO) 


The homogeneous boundary conditions (2.32)-(2.33) and the 
inversion of (2.48)-(2.50) give: 


ORIGINAU PAQf '® 
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4 R.-(a)p.^(mj2-va2) n»M^2 

^ f ^ r^A^(a) 

4 R4(a)P,-^ 

^ “ f ‘ ■'■ 

j2 I| R^(“)p^^ , , 

? rX2’'^mj2-A,‘!^){icp:':Tr ° 

4 

1 


= 0 


(2.52) 

(2.53) 

(2.54) 

(2.55) 


where 




qi(ot) = 


q2(“) 


q3(cx) 


G^(t)e^“^ dt 


(2.56) 




•v/c" 

G„(t)e^‘^^ dt 

(2.57) 

-v/E’ 


1'^ G3(t) e''"^ dt . 

(2.58) 

-/E" 



By solving the system of equations (2.51 )-(2.55) one can easily 
determine the unknown functions R-| (a) , . . . ,R^(ot) , A-j(a) in terms 
of G^(y), G 2 (y), G 3 (y)« The renaining mixed boundary conditions 
(2.34)-(2.36) give the following equations: 


lim, 

x-O'^ 


i f” 4 m.-x 

-00 



da = F,(y) 

t 


-/c < y < /c 


(2.59) 
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Itm - ^ 


---r — ipY e ^ dct 

.KPj’-lT 


a / V r o '♦ R4(a)P4^ni. 

Hrr n-v) ] a 2 a 2 

-00 C J I 

" x-lo'*' 2V I A^(a)(r.j2+a2)e ^ da « Fg(Y) 

->/c < y < \^ (2.60) 


x-^O"^ 2u 


r“ . R4(a)P4\ 


m.x 4,, 

— ™ p 1 q“ y« 
,tcPj-l) 


- lim 
x-K) 


c ^ X 

+ i K ^ a A^(a)e ^ e da = ^3(7) 


“v€* < y < /c 


(2.61) 


The Integrands In (2. 59) -(2. 61) are bounded for every value of 
a except when |a|-H» and y-^-t. Therefore one must extract the 
singular parts of the integrals VMhen |a|->«». The functions R-|(a), 
• ••» R/j,(a), A-j(a) may be written in the following form: 

Rj(a) *= i[Qj(a)qi(a) + Nj(a}q2(a) + Mj(a)q2(a)] 

(j « 1,...,4) (2.62) 

A-|(a) a B.j(a)o^(a) + Bg(a)q2(a) + B^(a)q2(a) (2.63) 

The functions Qj(a), Nj(a), Mj(a), (j=l,...,4) and B^(a), 
(i“l,2n3) are known from the solution of the system of equations 
(2.51) -(2. 55). Now using the asymptotic expansion of m. , 

J 

(j«l,...,4) and r^ for large values of |a| 
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Pn Pi 

mj(a) = -|a| (1 + 2^ “ 0^ +...)» (j=l >••• »4) (2.64) 

r-j(a) = -|a|(l + * (2.65) 


after some lengthy analysis, the singular parts of the inte- 
grals which appear in (2.59)-(2,61 ) can be separated and we 
obtain: 


‘/c G-j(t) 3 

' dt + E 


->/c 


t-y 


r/E 

J"”* -/E 


l<lj(y,t)Gj(t)dt 


= 2irF^(y) , -/c < y < }/E 

t n r^/E Gn(t) 3 

■' J-l 


( 2 . 66 ) 


= 2tt f F2(y) 


-/c < y < ^/E 


(2.67) 


^ Go(t) 3 

i ■ dt + E l<34(y,t)G.(t)dt 

0=1 




= 2 iTF 2 (y), < y < ^E . (2.68) 

The kernels k. .(y,t) which appear in equations (2.66)-(2.68) 

* J 

are given in Appendix B. 

To complete the formulation of the problem, one must also 
impose the single-valuedness conditions for the displacements 
and for the rotation of the normal, which are: 
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Gg(t)dt 

- 0 , 

(2.69) 










G3(t)dt 

= 0 , 

(2.70) 

-y/c 






X, 2 



[^^(t) + 

(-^) t G3(t)]dt = 0 . 

(2.71) 

iy/c 




The problem is thus reduced to tl<e solution of the singular 
integral equations (2.66)-(2.68) under the single-valuedness 
conditions (2.69)-(2.71 ) . 


a1 1 14* n-F ■fhfS liTf' 

*)/v/ivivr\/ii VI ullw iiiv 


cyi U I 


iqua 


4" n c 
~ 1/ I Ul lO 


Once the system of integral equations (2.66)-(2.71 ) is 
solved, then all field quantities can easily be computed in 
terms of the functions G.j(y), G2(y), G2(y). To do this, first 
the equations are normalized by defining: 


t = f/c r , -v^ < t < y^, “1 < X < 1 


(3.1) 


y = /c n j < y < /c , -1 < n < 1 , (3.2) 

x=»/c5,0<c,x<“’ (3.3) 


t) = H.(x) . (i = 1,2,3) 


(3.4) 


The functions 
let 

H^-(x) 


Hj(x), (i=l,2,3) are singular at x=+l. 


h^. (x) 

(x+1 )*^( 1-x) ^ 


(i=l,2,3) 


Therefore, 

(3.5) 
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where h^-Cr) are bounded and Holder-continuous In -1 < t < 1. 
Noting that the index of the problem k = (a+e) « 1 , and using 
the function-theoretic method described in [27] and [28], we 
have a ” 3 - ^ • Then, the integral equations can be solved 
numerically by using a Gaussian quadrature type formula [29-31], 


4. Asymptotic behavior around the crack tips. 


As stated earlier, one of the main objectives of this 
study is to show that by using a shell theory which adequately 
takes into account the effect of transverse shear strains, one 
can remove the discrepancy that exists between the classical 
shell theory and the elasticity solutions regarding the angular 
distribution of the stresses. Using the expressions given by 
(2.64) and (2.65), and the relation [28] 


1 

f 


-1 


h(t) i3x 


dx = (^ |^g| )^{h(l ) exp[i(3 - ^ jlj-)] 


+ h(-l)exp[-i(3 " f jly)] + (|3|-^) (4.1) 


around n=l and 5=0, the asymptotic expressions for the stress 
resultants are found to be: 


?hi(l) f“ oui 

^ ^ e“^'^l cos[3(l-n) - j]d3 (4.2) 

ZA'n J 




COS [3(l-n) - •r]d3 

/a 


(4.3) 
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N 


.h^d) ^3-315 


2v^ 

g h^O) 


e sin[e(l-n) - J ]d6 

/3 ^ 

u 

00 

[ /e sin[3(l-ri) - Jldp 


•^xx = hf^ ^ cosC3(l-n) 


2»^ 


M. 


= a 


,'y “ hx^ 


(l-v^) ^ f^^cos[8(l-n) - f]d6, 
v'^tT •' 

0 


M 


- i ilMH 

y.y ■ 12 a ^ J 


-elg 

^ sin[B(l-n) - jldp 


J_ il 

a 2v^ 


/3 sin[3(l-n) - f]d3, 


V.. = 




f3|€l 

/3 


sin[3(l-ri) " |-]d3 , 


ho(d f“ -eUl 

Vw = “=T ^“TT- cos[3(l-n) - f]d3 . 

^ >/27 J /3 

0 


If we now define the new coordinates. 


C = r sine , n-1 = r cose 


and use the following expression [32] 


(4.4) 

■]de 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 
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9 

(s > 0 » y > 0) * 

(4.11) 

(4.2)-(4.9) become; 


h-j(l) 1 ^ 0 i.ie.Jn 

N = - J_ [- 7 sin ^ + 4 sin , 

(4.12) 

^ 1 ^^^ r 7 • 0 i 50n 

(4.. 13) 

„ ; . [|cos| + {cosf ] , 

(4.14) 

u hn(l) T A T *^6 - 

f A n f* \ 

(4f. 10 ; 

i. h«(l) _ 7 . 0 1 • 50 T 

M s - -.y — = [- -r sin 7 " j sin J » 

yy 12= 2 y 2 r ^ ^ 

(4.16) 

h ^p(i) _ 3 6 1 50 n 

M = - t 5 T [ f COS ^ 4 cos -^ ] , 

xy 12a 2 ,^ 4 2 a 

(4.17) 

- 1 ^ 3 ( 1 ) 0 

X ^ 2 

(4.18) 

' 113 ( 1 ) g 

V = — sin %■ . 

(4.19) 


It should first be noted that the new coordinates defined by 
( 4 . 10 ) are not polar coordinates for ’'specially" orthotropic 
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materials. For isotropic materials 1t is seen that the dis- 
tributions given by (4,12-4.19) are identical to those obtained 
from the plane stress and anti-plane elasticity solutions [18], 
[33]. 


5. Results and Discussion 

After solving the system of equations (2.66)-(2.71) , the 
discrete values of the unknowns G-j(t), GgCt), G^Ct) can be used 
to determine any desired field quantity in the shell. The 
mode II and mode III stress intensity factors are defined as: 

k2(X3) = lim /2Txpa7 
X2->« 

kJ^) = lim 0,3 (0, X«, X3) (5,2) 

C 10 CO 

Using expressions (4.12)-(4.19) and those given in Appendix A 
with 6=0, the stress intensity factors can easily be expressed 
in terms of the end values of the unknown functions. In the 
examples given, three types of loading are considered. 

a) In-plane uniform shear loading: 

If denotes the magnitude of uniformly distributed 
shear stresses through the thickness, the external loads can 
be expressed as: 

I 

Ni 2(0, Xg) = -N^2 "" ’ 

^12^^* ^2^ ” 0 , 

V^(0, X2) = 0, -a < X2 < a (5.3) 
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using the dimensionless quantities defined in Appendix A, the 
functions F^, Fj, Fj that appear in (2.66)-(2.68) are foun 

to be: 

F,(v€ n) = - ^ • n) = 0 , FjC.^ n) = 0, 

-1 < n < 1 . 

In this case the stress intensity factors will be normalized 
with cf—v^ • Thus* we have. 


'm 


(< 0 ( 0 ) h-|(l) £ 


a„>/a 

m 


2 ’ 


(5.5) 


— ^ S “ n' _ O-t * 


'tm 


m 


2 2a 
m 


(5.6) 


= h,(l) . 


S'" O 

m 


2 a, 


m 


(5.7) 


b) Uniform Twisting Moment 

Let Of denote the maximum of linearly distributed shear 
stresses through the thickness. Then, the external loads are 

^12^^’ ^2^ ” ^ * 

h^ 

M - j 2 ( 0 » ^ 2 ^ ~ “^12 " " 6 *^1 ’ 

V^(0, Xg) = 0 , -a < X 2 < a , 


or 
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(5.9) 


F-j(v^ n) ® 0, ^2^"^ n) “ - n) =* 0 , 


-1 < n < 1 


Normalizing the stress intensity factors with respect to , 
we obtain: 


mt 


_ kgCO) _ h^(l) g 


at»^ 


2 0j 


(5.10) 


'tt 


‘st 


_ kg(h/2)-k3(0) hjd) £ ^ 


^o(O) O D 


at 2a 


(5.11) 


(5.12) 


c) Uniform Transverse Shear Loading: 

If 0 g denotes the maximum of the parabolical ly distributed 
transverse shear stresses, the external loads can be expressed 
as: 

\ 

N,2(0, Xj) = o’, 

Xj) = 0 , 

y^(0, ^2^ ~ ’ ’■9*'X2<s (5,13) 

or 

F^Cv/^n) = 0. F 2 (v/c n) = 0, n) = - |x“ » 

/c 

‘1 < n < 1 . (5.14) 
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In this case, normalizing the stress intensity factors with 
respect to , we obtain: 


k ® 
ms 


B 

“ts 


l<2(0) h,(l) g 


OjA" 


^ °s ’ 


k3(h/2)-kg(0) _ hgd) E I, 


a 

s 


2 cTg 2a * 


(5.15) 


(5.16) 


ag/a s 


(5.17) 


The stress intensity factor ratios are calculated for different 
shell geometries. The following examples are considered; 


a) 

b) 

c) 

d) 


a cylindrical shell with an axial crack, 
a circumferentially cracked cylindrical shell, 
a spherical shell with a meridional crack, and 
a toroidal shell containing a crack at different locations, 


In all these examples, the material is assumed to be isotropic, 
with v-0.3. To give also some idea about the effect of ortho- 
tropy on the results, an example is considered for an axially 
cracked cylinder. The elastic properties of titanium used in 
these calculations are as follows: 

= 1.039 X 10^Vm2 , 

E 2 = 1.434 X 10^^N/m2 , 

v-j = 0.1966, 

Vg =* 0.2714, 


no 


G^2 “ 4-675 X 10^° N/m2 , 

“ 4.956 X 10^° N/in2 , 

aVc 


The results are given in tables 1-34. Some of the results are 
also displayed in figures 3 and 4. First, it muj^t be pointed 
out that for toroidal shells, according to the crack location 
or are taken as negative to describe negative curva- 
tures. Secondly, as x-j or Xg approach zero, the results given 
in the tables approach those obtained for flat plates £18J. 
From tables 1-34 and figures 3-4 the following trends can be 
observed: 

- If the crack surfaces are loaded with in-plane shearing 
forces, the coupling stress intensity factor ratios k^^ and 
kgj^ are small compared to As Xg increases, i.e. as the 


"I 

1 nw I caaca ^ 


mm 


1 1 ^ 
( I aw 


■f irt na 
I ii^t cuavsa • 


n« 

VI I 


VI iw 


rk4*hov< hanri 

wi 1 wi I ivft I IW 9 


a/h which expresses the thickness effect, does not affect the 
results significantly. 


- Under uniform twisting moment applied to the crack sur- 
faces, the stress intensity factor ratio k^.^ becomes dominant, 
and the other two components k^^^ and k^^ are small in compari- 
son. In this case, however, k^^ does not vary significantly with 
X 2 » and takes smaller values as a/h increases. 

- For uniform transverse shear loading, the stress intensity 

factor ratio k^^g is very small, however, k^g is no longer 

negligible. For this case also, the dominant stress intensity 

factor ratio is k^^. Tne stress intensity factor ratio k^_ 
ss ss 

does not vary significantly with curvature however, and decreases 
as the thickness of the shell increases. 
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- Figure 3 shows the comparison of for three shell 
geometries, assumes its smallest values for a cylinder 
with an axial crack, and assumes larger values for a circum- 
ferentially cracked cylindrical shell. The results for a 
spherical shell are even larger. 

- The results given in tables 28-31 for toroidal shells show 
the same trends for the stress intensity factors as the cylin- 
drical shells. It must be pointed out that, for this case, the 
component k^^ remains almost unaffected of the location of the 
crack. 

- The effect of orthotropy on the results is given in 
tables 33-34 and figure 4. As it may be observed, k^^ is the 
component that is affected most by material orthotropy. If thv» 
axes of orthotropy are subjected to a 90° rotation, the results 
change significantly. For example, for E-j/Eg < 1 are 
smaller than those obtained for isotropic materials, whereas 
for E-j/Eg > 1» are larger than the isotropic results. 
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Appendix A 

Dimensionless quantities used in th... derivation: 


^ 1 .H 
" ^ 


, y = /c — , z 


.'3 


/- ^1 1 ^2 W 

v=--^.w = - 


= v/c 61 , 3 y = -^ ^2 * " a^ 


'^xx "" <^11/^^ ’ ^yy " ^ ^22/^ ’ ‘^xy "" ^^12/^ 

“ ) 0 ^,^ “ /c" ^ 23 /^ 


^xz '^IS 


yz 


N,x = N„/chE , Nyy = cN22/hE , N^y = N^j/hE 


”xx = M,,/ch2E . Myy = cMgj/h^E , 


= V-j//c h B , Vy = /c V 2 /I 1 B 
^ 1 “ = 12(l-v2) ’ >- 2 “ " 12(1-''^) 


^12 = 12(l-v2) , x“* = 12{l-v2) |r . ^ = gf^ 


(A.l) 


(A.2) 


(A. 3) 

(A. 4) 
(A. 5) 
(A.6) 
(A. 7) 
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Expressions 
ki2(y>i') ~ 

k-igCy.t) « 
k2i(y.t) = 

k23(y.t) = 
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Appendix B 

of the kernels which appear in the integral equations; 
•j [2ia E nijQj(a)+l]sina(t-y)da (B.l) 


li 


H 

2ia E m.N.(a)sina(t-y)da 

1 j j 

” 4 

a E nijMjCa)cosa(t-y)da 


C (a2+r^2)B^(„)]s1no(t-y)dci 

fr„, 1-. 5 


(B.2) 


CB.3) 


r ri-v i "i i '".i 

>■1^“ ? ’> 


(B.4) 


i - (a^+ri^)B 2 (ci) - si no (t-y)dg (B.5) 

M^(a)p,.2m.^ 

■J - fej 


(a2+r-j^)B2(a)]coSa(t-y)da 


CB,6} 
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kg-lCyjt) - 

-i 


i<33(y»t) = 
+ 


’fh- ^2 K E Qi(»)Pi^'"i 

■J. ' 




K il^yl C( B^Ca)]cosaCt-y)da 

f” 4 N.Ca)Pi^ni. 


CB.7! 


1 K a B 2 (a)]cosa(t-y)da 

f 4 M.(a)p.3m. 

- j *= I (X22m%x!,i2^(Kp:-TT 


CB.8) 


< B^Ca) - l]sina(t-y)dot 


C.B.9) 
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Table 1. Stress Intensity factor ratios kipm for an iso- 
tropic cylinder containing an axial crack, v=0.3. 
(Uniform in-plane shearing) 


^1 

mmsm 



a / h =4 

0.0 

1.000 

1.000 

1.000 

1.000 

0.25 

1.003 

1.003 

1,003 

1.003 

0.50 

1.012 

1.011 

1.011 

1.011 

0.75 

1.026 

1.024 

1.024 

1.024 

1.0 


1.040 

1.039 

1.039 

1.5 


1.078 

1.076 

1.075 

2.0 



1.114 


3.0 






Table 2. Stress intensity factor ratios k^f^ for an iso- 
tropic cylinder containing an axial crack, v=0.3. 
(Uniform in-plane shearing) 


M 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

- 0.005 

- 0.002 

- 0.001 

- 0.001 

0.50 

- 0.012 

- 0.005 

- 0.002 

0.000 

0.75 


||||Q|gg||||||H 

- 0.001 

0.003 

1.0 

bhh 

- 0.006 

0.003 

0.008 

1.5 


0.002 

0.017 

0.027 

2.0 



0.037 

0.051 

3.0 




0.103 
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Table 3. Stress intensity factor ratios ksm for an isotropic 
cylinder containing an axial crack, v=0.3. (Uniform 
in-plane shearing) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

0.005 

0.004 

0.003 

0,003 

0.50 

0.020 

0.014 

0.012 

0.011 

jj^glll 

0.042 

0.030 

0.026 

0.024 

1.0 


0.051 

0.043 

0.039 



0.101 

0.086 

0.078 

2.0 



0.137 


3.0 




0.222 


Table 4. 


Otrtibb 


In tens! i.y 
cylinder containing 
twisting moment) 


I aCuui 


an 


ratios Knrc lor an 
axial crack, v=0.3. 


1 sotropi c 
(Uniform 


^1 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 


- 0.001 

- 0.001 

0.000 

0.50 

HHHI 

- 0.002 

- 0.002 

- 0.001 

0.75 

- 0.007 

- 0.004 

- 0.003 

- 0.002 

1.0 


- 0.005 

- 0.004 

- 0.003 

tsm 


- 0.007 

- 0.005 

- 0.004 

2.0 



- 0.006 


3.0 




- 0.005 
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Table 5. Stress intensity factor ratios for an 

isotropic cylinder containing 'n axial crack, 
v=0.3. (Uniform twisting moment) 



a/h=l 

a/h=2 

a/h=3 

a/h®4 

0.0 

0.523 

0.354 

0.274 

0,228 

0,25 

0.523 

0.354 

0.274 

0.228 

0.50 

0.522 

0.353 

0.274 

0.227 

0.75 

0.521 

0.353 

0.274 

0.227 

1.0 


0.353 

0.274 

0.227 

mm 


0.352 

0.273 

0.227 

2.0 



0.273 

0.227 

3.0 




6.226 


Table 6. Stress intensity factor ratios kst for an isotropic 
cylinder containing an axial crack, v=0.3. (Uni- 
form twisting moment) 


^1 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 

-0.070 

-0.092 

-0,094 

-0.091 

0.25 

-0.070 

-0.092 

-0.094 

-0.091 

0.50 

-0.070 

-0.091 

-0.094 

-0.091 

0.75 

-0.069 

-0.091 

-0.094 

-0.091 

1.0 


-0.091 

-0.093 

-0.091 

1.5 


-0.089 

-0.092 

-0.090 

2.0 



-0.091 

-0.089 

3.0 




-0.087 
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Table 7. Stress Intensity factor ratios kjps for an iso- 
tropic cylinder containing an axial crack, v=0.3 
(Uniform transverse shear loading) 



a / h®l 

a / h «2 

a / h «3 

a / h -4 


0.000 

0.000 

0.000 

0.000 


n^m^gniii 

“ 0.005 

“ 0.006 

“ 0.007 

0.50 

“ 0.021 

“ 0.021 

- 0.024 

“ 0.028 

0.75 

“ 0.045 


“ 0.052 

“ 0.061 

1.0 


“ 0.076 

“ 0.088 

“ 0.103 



“ 0.149 

“ 0.174 

“ 0.204 

?,.o 



“ 0.264 

- 0.310 

3.0 




“ 0.495 


Table 8. Stress intensity factor ratios k^s an iso- 
tropic cylinder containing an axial crack, v=0,3 
(Uniform transverse shear loading) 


^1 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.47 

1.26 

2.12 

3.03 

0.25 

0.47 

1.26 

2.12 

3.03 

0.50 

0.46 

1.26 

2.12 

3.03 

0.75 

0.46 

1.25 

2.11 

3.02 

1.0 


1.24 

2.10 

3.00 

1.5 


1.19 

2.04 

2.93 




1.95 

2.81 

3.0 




2.52 


- 34 - 
































































ORIGlMAt 

OF POOR QUAUfi 


Table 9. Stress intensity factor ratios kss for an isotropic 
cylinder containing an axial crack, veQ.3. (Uni- 
form transverse shear loading) 


^1 

a/hrl 

a / h *2 

a / h «3 

a / h -4 

0.0 

1.68 

2.34 

2.98 

3,62 

0.25 

1.68 

2,34 

2.98 

3.62 

0.50 

1.67 

2.34 

2.98 

3.62 

0.75 

1.66 

2.33 

2.98 

3.61 

1,.0 


2.32 

2.96 

3.60 

1.5 


2.28 

2.92 

3.54 

2.0 



2.84 

3,46 

3.0 




3.24 


Table 10. Stress intensity factor ratios for a cir- 
cumferentially cracked cylindrical shell* 


\2 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

1.000 

1.000 

1.000 

1,000 

0.25 

1.003 

1.003 

1.003 

1.003 

0.50 

1.013 

1.012 

1.012 

1.012 

0.75 

1.031 

1.027 

1.026 

1.026 

1.0 


1.047 

1.045 

1.045 

1.5 


1.104 

1.097 

1.095 

2.0 



1.165 

1.159 

3.0 




1.320 
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Table 11. Stress intensity factor ratios k^m for a cir- 
cumferentially cracked cylindrical shell. 


\z 

a / h=l 

a / h =2 

a / h ®3 

a / h =4 

0.0 

0.000 

0,000 

0.000 


0.25 

-0.005 

nmggnii 

0.000 

0.001 

0.50 

■mgmuii 


0.002 

0.006 

0.75 


-0.002 

0.010 

0.017 

1.0 


0.004 

0.023 

0.035 

ingm 


0.030 

0.068 

0.092 

2.0 



0.139 

0.178 

3.0 




0.437 


Table 12. Stress intensity factor ratios kgiji for a cir- 
cumferentially cracked cylindrical shell. 


Xz 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

0.014 

0.010 

0.008 

0.007 

0.50 

0.056 

0.038 

0.032 

0.028 

0.75 

0.125 

0.034 

0.070 

0.062 

1.0 


0.148 

0.123 

0.109 



0.328 

0.272 

0.242 

2.0 



0.479 

0.426 

3.0 




0.949 
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Table 13. Stress intensity factor ratios !<„,+ for a 

circumferentially cracked cylindrical shell. 
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Table 15. Stress intensity factor ratios kst for a cir- 
cumferentially cracked cylindrical shell. 



Table 16, Stress intensity factor ratios kjy,.. for a cir^ 
cumferential ly cracked cylindrical shell. 
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Table 19. Stress Intensity factor ratios for a spheri- 
cal shell with a meridional crack. 


0.0 


0.25 


0.50 


0.75 


a / h=l 


1.000 


1.006 


1.026 


1.058 


a / h=2 


1.000 


1.006 


1.024 


1.051 


1.086 


1.172 


a / h=3 


1.000 


1.006 


1.023 


1.049 


1.082 


1.162 


1.252 


a / h-4 


1.000 


1.006 


1.023 


1.048 


1.081 


1.158 


1.244 


1.418 


Table 20. Stress intensity factor ratios ktm i"or a spheri 
cal shell with a meridional crack. 


0.000 


- 0.009 


- 0.021 


- 0.030 


a / h=2 

a / h=3 

0.000 

0.000 

- 0.004 

- 0.002 

- 0.005 

0.003 

- 0.001 

0.016 

0.011 

0.037 

0.056 

0.105 


0.205 


- 0.001 


0.008 


0.025 


0.052 


0.136 


0.254 


0.581 


0 - 
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Table 21. Stress intensity factor ratios ksm for a spheri- 
cal shell with a meridional crack. 


X 2 

a/h=l 

a/h=2 

a/h-3 

a/h=4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

0.014 

0.010 

0.009 


0.50 

0.077 

0.053 

0.045 

0.040' 

0.75 

0.169 

0.115 

0.096 


1.0 


0.200 

0.166 

0.149 

1.5 


0.436 

0.361 

0.322 

2,0 



0.626 

0.558 



3.0 






Table 22. Stress intensity factor ratios k^t for a spheri- 
cal shell with a meridional crack. 


X 2 

a/h=l 



a/h==2 

a/h=3 

a/li-4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

-0.003 

-0.002 

-0.001 

-G.OOl 


-0.008 

-0.004 

-0.003 

-0.002 

0.75 

-0.012 

-0.007 

-0.005 

-0.004 

1.0 


-0.010 


-0.005 

Dgjll 


imimiii 


iimgsg 

2.0 



imiggiiii 

-0.009 

3.0 




-0.012 
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Table 23. Stress intensity factor ratios ktt for a spheri 
cal shell with a meridional crack. 


a/h-1 

0.523 

0,523 

0.521“ 

0.519 


a/h=2 

0.354“ 

“oT35^ 

0.353 

0.352“ 

■^OsT 

0.350 


a/h=3 

0.274 

0.274 

0.274“ 

0.27T 

0.273 

0.272“ 

0.271 


a/h=4 

0.22i 

0.227 

0.227 

0.227 

0.227 

0.226 

0.226 

0.226 


Table 24. Stress intensity factor ratios kst for a spheri 
cal shell with a meridional crack. 


a/h=2 


a/h=3 




DJIil 


-0.091 


-0.09 
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Table 25. Stress intensity factor ratios k^s a spheri- 
cal shell with a meridional crack. 


m 

a / h=l 

a / h =2 

a / h =3 

alhr ^ 

0.0 

0.000 

0.000 

0.000 

0,000 

0.25 



- 0.025 

- 0.029 

0.50 

- 0.085 


- 0.094 

- 0.108 

0.75 



- 0.199 

- 0.230 

1.0 


- 0.282 

- 0.325 

- 0.379 

1.5 



- 0.563 

- 0.669 

2.0 



- 0.691 

- 0.844 

3.0 




- 0.843 


Table 26. Stress intensity factor ratios k-ts for a spheri- 
cal shell with a meridional crack. 


^2 

a / h=l 

a / h =2 

a / h =3 

a / h -4 

0.0 

0.47 

1.26 

2.12 

3,03 

0.25 

0.46 

1.26 

2.12 

3,03 

0.50 

0.45 

1.25 

2. 11 

3.00 

0.75 

0.42 

1.20 

2.05 

2.94 

1.0 


1.12 

1.95 

2.81 

1.5 


0.86 

1.58 

2 . "3 5 

2.0 



1.15 

1.76 

3.0 




0.84 
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Table 27. Stress intensity factor ratios kgs for a spheri- 
cal shell with a meridional crack. 



Table 28. The stress intensity factor ratios for the crack 
configuration shown in figure (2. a). Toroidal 
isotropic shell, with v=0.3. 
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Table 29. The stress intensity factor ratios for the crack 
configuration shown in figure (2.b). Toroidal 
isotropic shell, with v=0.3. 




Table 30. The stress intensity factor ratios for the crack 
configuration shown in figure (2.c). Toroidal 
isotropic shell, with v=0.3. 
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Table 31. The stress intensity factor ratios for the crack 
configuration shown in figure (2.d). Toroidal 
isotropic shell, with v=0.3. 


a/h=l 

Xl=0.332 

X2=0.575 


0.073 


0.45 


1.64 


aAi=2 

Xl=0.664 

A2-1.150 



0.06 


0,352 


-0,091 


0.247 


1.11 


2,08 


a/h=3 
X 1=0. 996 
X2-1.725 


1,100 


.08 


-0,312 


0.007 


0.273 


93 


0.485 


1.50 


2.18 


a/h=4 

Xi=l,328 

X2«2.300 


1.164 


-0,488 


0,008 


0,227 


0,700 


1,57 


1.99 


Table 32. The stress intensity factor ratios for the crack 
configuration shown in figure (2.e). Toroidal 
isotropic shell, with v=0.3. 


Ri/R2= 1/5 
X2=0.514 

i 

Ri/R2= 1/12 
X2=0.332 

Ri/R2= 1/20 
X2=0.257 

1.063 

1.056 

1.054 


kabuk 
Ri/R2= 0 
X2=0 


1.051 


0.081 


-0.006 


0.352 


s 


1.23 


1.23 


















































Table 33. Stress intensity factor ratios for an orthotropic 
(titanium) axially cracked cylindrical sh ell . 
(a/h=3, E3^/E2=0.725, Xq= [I2(l-v2)]l/'*a/*^^) 



























































Table 34. Stress intensity factor ratios for an orthotropic 
(titanium) axially cracked cylindrical she ll . 
(a/h=3, E3^/E2=1.380, Xq= [12(1-v^)] 1A a/i^Jh) 
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Companson of the stress intensity factor ratios 
Kim for cylindrical and spherical shells 




